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Abstract 


The comfort in car riding largely depends on its suspension system. Isolation of passenger 
compartment from the vibrations due to road undulations is the primary purpose of the suspension 
system. Two-dimensional as well as three-dimensional dynamic models of the suspension system 
have been developed with front and rear suspensions as of independent type. Springs and 
dampers of the suspension system can be chosen to make the isolation maximum, while satisfying 
other important criteria of the performance of the system. The optimization for the design of 
suspension system is carried out using genetic algorithm — a optimization technique based on the 
principles of natural genetics. Various idealized road conditions are used with different objectives 
related to different motions of the sprung mass of the car. to get the optimal solution. Vertical 
jerk experienced by the passenger is kept below a certain level by using a constraint on its value. 
Natural frequencies of the subsystem of the suspension system imposes other constraints on the 
optimal solution. A realistic road is modeled using a polyharmonic function. ISO standards are 
used to find the level of comfort based on the acceleration levels on different frequencies. The 
solution found in all the cases are better than the present design of TELCO, Pune. Success of 
genetic algorithm in finding the near optimal solution in highly complex, multimodal and 
disjointed feasible search space of the optimal suspension design problem suggests their 
immediate application in other problems of engineering. 
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Chapter 1 


Introduction 


In today’s competitive world the need of the time is to use the available resources in a most 
effective manner. This means that, whenever possible, the maximization of the performance of any 
system should be achieved. In engineering design problems, the measure of the performance 
varies from problem to problem. The maximization of performance may be in terms of reducing 
overall cost of the production or it may be in terms of maximization of the profit for a product 
|7]. This should be achieved with other related aspects of the problem in mind. In simpler words, 
the performance should be optimal taking care of all conflicting aspects of the problem. 
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The modern-day vehicle is a complex assembly created to meet a variety of design 
objectives, including fuel economy, handling, performance, safety, passenger comfort, 
driveability, and ride quality. No single vehicle can meet all of these objectives maximally. As a 
result, designers make compromises on some objectives while emphasizing some other design 
objectives. Suspension system of any vehicle has to meet various requirements for the proper and 
smooth functioning of the vehicle. Suspension is supposed to support the vehicle weight and to 
provide proper handling characteristics of the vehicle. Most importantly, it attempts to isolate the 
passenger compartment of the vehicle form vibrations due to road roughness. The effect of road 
vibrations experienced by the passenger interferes with the comfort during the ride. The 
perception of the ride comfort is a subjective matter but in general vehicle designed for good 
comfort should have a suspension system which reduces the transmissibility of the road 
undulations to its occupants. 

In this thesis, we attempt to find an optimal design of suspension system of a car with the comfort 
of the passenger as the objective. The dynamic model of the suspension system is achieved 
mathematically to simulate the problem. Firstly, a two-dimensional model is developed and used 
for solving the problem. Thereafter, a three-dimensional model is developed for more detailed 
study of a car’s motion and for getting a better design. Different motions of the car are taken into 
consideration and objective functions are formulated. The jerk experienced by the occupants is 
taken as a constraint. Natural frequencies of the front and rear subsystems of the suspension 
system also impose some constraints to the optimal problem. 

To solve the optimization problem, the genetic algorithm is used. Genetic algorithms are new yet 
potential techniques based on the principles of natural genetics [6,12,14]. Population-based 
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approach of the genetic algorithms enables them to find the global optimum solution in most 
complex problems. Operators of these algorithms have remarkable ability of combining 
information of two partial solutions to get a better solution [14]. In this thesis, GA technique is 
tested well for finding optimal solution for different objectives and for different road excitations. 
ISO provides standards on human exposure to whole-body vibrations [15,20]. Some limits have 
been discussed in these standards for the comfort in such conditions. Motion of a car over a rough 
road gives this kind of exposure to its occupants. An attempt has been directed towards getting 
an optimal design in accordance of these standards. 

In Chapter 2, we develop the theory of the suspension design procedure, with two as well as 
three-dimensional conditions. Chapter 3 gives an overview of the genetic algorithms discussing its 
working principle and properties of its basic operators. Some previous applications of GAs to 
engineering problems are also presented. The optimal design problem formulations for different 
objectives related to the suspension design are presented in Chapter 4. This chapter also presents 
the formulation of the constraints related to natural frequencies of subsystems of the suspension 
system. Chapter 5 presents the results of different simulation runs of a traditional optimization 
method and of a GA for finding optimal solutions. Simulation of the car motion over a realistic 
road and subsequent analysis in accordance with ISO 2631 are also discussed. Conclusions and 
extensions of the present work are presented in Chapter 6. 
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Suspension Systems: 
Theory And Modeling 


The suspension systems, like other parts of the vehicle, are designed to meet specific 
requirements. These requirements consider the number of passenger to be carried, the useful load, 
the type of the ride desired and the required handling characteristics. The primary function of the 
suspension system is to isolate the structure, as far as possible, from shock loading and vibrations 
due to irregularities of the road surface. Secondly, it must do this without impairing the stability, 
steering or general handling qualities of the vehicle. The primary requirement is met by the use of 
flexible elements and dampers while the second is achieved by use of mechanical linkages [17]. As 
is generally true, when one suspension property is improved, the quality of another property is 
reduced. In this thesis, we concentrate on designing the former for given road conditions. 
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2.1 Vehicle Dynamics 

Progress on understanding the ride motion of a vehicle is made by applying the principles 
of vibration theory. However, putting this theory into practice is not so straightforward, since it 
includes several suspension design issues. The designer has to achieve good ride comfort for the 
driver and passengers, acceptable control of the body attitude and adequate control of the 
dynamic tyre loads. The diameter of the tyre, size of the contact patch between tyre and the road, 
the rate of tyre acting as a spring, and weight of the wheel and axle assembly affect the 
magnitude of the shock transmitted to the axle, while amplitude of wheel motion is influenced by 
all these factors plus the rate of suspension springs, damping effect of the shock absorbers and 
weight of sprung mass and unsprung mass. The unsprung mass can be loosely defined as that 
between the road and the main suspension springs, while the sprung mass is that supported on 
these suspension springs [10,17]. 

Dampers in suspension system have a two-fold function. First, they are for reducing the 
tendency for carriage unit to continue to bounce up and down on its springs after the disturbance 
that causes the initial motion has ceased. Secondly, they prevent excessive building-up of the 
amplitude of the bounce as a result of periodic excitation at a frequency identical to the natural 
frequency of the vibration of the system. Thus, for maximum ride comfort, proper spring stiffness 
and damper characteristics must be chosen. 
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2.1.1 Ride Characteristics 

The vehicle is a dynamic system but only exhibits action response to an excitation input. The 
response properties determine the magnitude and direction of vibrations imposed to the 
passenger compartment and ultimately determines the passenger’s perception of the vehicle. The 
understanding of ride involves the study of following [11]: 

• Ride excitation sources 

• Basic mechanics of the vehicle vibration response 

• Human perception and tolerance of vibration 



2.1.1.1 Excitation Sources 

There are multiple sources from which vehicle ride vibrations may be caused. These are 
generally divided into factors such as road roughness and on-board sources. The road roughness 
comes from the roughness of the road, while the on board sources arise from reciprocating 
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components of the car. However, m this thesis we shall concentrate on the factors related to road 
roughness only. 

Road roughness encompasses everything from potholes to the ever present random 
deviations reflecting the practical limits of precision to which the road surfaces can be constructed 
and maintained. Roughness is described by the elevation profile along the wheel tracks over which 
the vehicle passes. Road profiles fit the general category of the “broad-band random signals” and 
can be described by the profile itself or its statistical properties [11]. 

Like any random signal, the elevation profile measured over a length of road can 
be decomposed using the Fourier Transform into a series of sine waves varying in their 
amplitudes and phase relationship. In chapter 5, we simulate the road by assuming it as an 
addition of sine waves defined in varying wavelengths. 

2.1.1.2 Vehicle Response Properties 

A systematic treatment of the vehicle as a dynamic system begins by identifying 
basic supporting structure of a vehicle on its suspension system, i.e., the motion of the body and 
axles. At low frequencies the body, which is considered to be the sprung mass portion of the 
vehicle, moves as an integral unit of the suspension. This is like rigid body motion. The axles and 
associated hardware, which form the unsprung masses also moves, as a rigid body and 
consequently impose excitation forces on the sprung mass. 

The dynamic behavior of the vehicle can be characterized most meaningfully by 
considering the input-output relationships. This input may be any one of the excitations discussed 
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earlier. The output most commonly of interest will be the vibrations on the body. The ratio of 
output and input amplitudes represents “gain” for the dynamic system. The term transmissibility 
is often used to denote this gain. Transmissibility is a non-dimensional ratio of the maximum 
amplitude of vibration of the sprung mass to the amplitude of excitation for a system in forced 
vibration. Although this ratio may be calculated for amplitudes of forces, displacements, 
velocities or accelerations, we have considered amplitude of linear and angular displacements. 

2.1.1.3 Perception of Ride 

Theoretical assessment of the ride vibrations must deal with the issue of how a 
ride is perceived. For that purpose, one must first attempt to define ride. Ride is a subjective 
perception, normally associated with the level of the comfort experienced when traveling in a 
vehicle. Therefore, in its broadest sense, the perceived ride is the cumulative product of many 
factors. The tactile vibrations transmitted to the passenger’s body through seat, and at the hand 
and the feet, are factors most commonly associated with the ride. Yet it is often difficult to 
separate the influences of acoustic vibrations (noise) in the perception of a ride. Additionally, the 
general comfort level can be influenced by seat design and its fit to the passenger, temperature, 
ventilation, interior space, hand holds, and many other factors. These factors may contribute to 
what might be termed as the “ride quality” of a vehicle [11]. Some of the above factors, such as 
vibrations, can be measured objectively, while others, such as seat comfort, are still heavily 
dependent on the subjective evaluation methods. Measurement of vibrations in terms of 
transmissibility of the vertical and angular displacement is taken as criterion for the assessment of 
the ride comfort in this thesis. 
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2.2 Type Of Suspension Systems 

Suspension systems have been changed and refined as passenger automobile has 
developed. Design objectives differ between luxury sedans, performance vehicles, small compact 
vehicles, and light trucks. Design changes for different engine positions, engine mountings and 
drive. Here, we discuss different suspension system commonly used in passenger cars. 

2.2.1 Front Suspension 

Front suspension designs have been developed from relatively rugged solid-axle 
designs to the modem lightweight, high-strength, stmt-type independent designs. These have been 
upgraded with added linkage. 

Dependent-type, front-suspension systems with solid axles were used on early 
domestic vehicles [10]. On a dependent suspension the motion at one wheel causes movement at 
other wheel. The solid-axle front suspension is used on vehicles that go over rough roads at low 
speeds. It is the strongest type of the front suspension system (Figure 2. 1). 

The independent front suspension system grew from the demand for improved nde 
quality and improved handling at increased speeds. In independent suspensions, motion of one 
wheel does not cause movement at the other wheel. The independent front suspension primarily 
used the long-short arm design (Figure 2.2). Independent front suspension designs on some 
vehicles took the form of either a leading or a trailing arm design. Today, an increased number of 
vehicles uses a strut-type of design (Figure 2.3). The stmt design provides good ride, good 
handling, and increased underhood space. 
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Figure 2.1: Typical Dependent Front Suspension (taken from [10]) 



Figure 2.2: Long-short Arm Type Independent Front Suspension (taken from [10]) 


Figure 2.3: Strut Type Front Suspension (taken from [10]) 
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We will consider independent type of front suspension for modeling the suspension system , as 
this type of suspension is applied in most of the cars designed for comfort ride. 

2.2.2 Rear Suspension 

Most front-engine, rear-wheel-drive vehicles use a simple dependent rear 
suspension. Rear-wheel-drive independent suspension is much more complex and expensive. 

The move to front-engine, front-wheel-drive has allowed more flexibility in rear 
suspension design. By moving the drive train to the front, only ride control and breaking reactions 
are to be controlled by the rear suspension. This has led to the use of simplified dependent 
suspensions (Figure 2.4), semi-independent suspension (Figure 2.5), and independent rear 
suspension (Figure 2.6) . 

Independent type of rear suspension system is used in next sections for the modeling of the 
suspension system. 

2.3 Suspension Modeling 

At the most basic level, all vehicles share the ride isolation properties common to 
a sprung mass supported by primary suspension systems at each wheel. The essential dynamics 
can be represented by a quarter car model (Figure 2.7). It consists of a sprung mass supported on 
a primary suspension, which is in turn connected to the unsprung mass of the axle. The suspension 
has the stiffness and damping properties. The tyre is represented as a simple spring. 
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iRoad Profile 



Figure 2.7 Quarter Car Model 
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2.3.1 Two-dimensional Model 

The simple mechanics of quarter car model do not fully represent the rigid body 
motions that may occur on a car. Because of the longitudinal distance between axles, it is a multi- 
input system that responds with pitching motion as well as with a vertical bounce. Depending on 
the road and speed conditions, one or the other motions may be largely absent. A two- 
dimensional dynamic response is essential for the analysis of this motion. For this, the suspension 
system can be modeled as shown in Figure 2.8. It consists of a sprung mass supported on the two 
axles by means of primary suspension, that is, a suspension coil spring and a shock absorber 
(damper). Each axle contains a unsprung mass that is supported only by the tyres. 

To define the model completely the following parameters are to be specified: 

Sprung mass (M s ) 

Front and rear unsprung mass( M fu> , M m ) 

Front and rear coil stiffness (K fs , , K rs ) 
front and rear damper coefficient ( a/,, a r ) 

Front and rear tyre stiffness (Kp , K n ) 

Axle to axle distance (L) 

Polar moment of inertia of the sprung mass (7) 

The two-dimensional model of the car suspension system has four degrees of freedom, namely 

• Vertical motion of front unsprung mass( q\) 

• Vertical motion of sprung mass(^) 

• Angular motion of sprung mass(^ 3 ) 

• Vertical motion of rear unsprung massO^) 
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2.3.1 Two-dimensional Model 

The simple mechanics of quarter car model do not fully represent the ngid body 
motions that may occur on a car. Because of the longitudinal distance between axles, it is a multi- 
input system that responds with pitching motion as well as with a vertical bounce. Depending on 
the road and speed conditions, one or the other motions may be largely absent. A two- 
dimensional dynamic response is essential for the analysis of this motion. For this, the suspension 
system can be modeled as shown in Figure 2.8. It consists of a sprung mass supported on the two 
axles by means of primary suspension, that is, a suspension coil spring and a shock absorber 
(damper). Each axle contains a unsprung mass that is supported only by the tyres. 

To define the model completely the following parameters are to be specified: 

Sprung mass (M s ) 

Front and rear unsprung mass( M fu> , M m ) 

Front and rear coil stiffness (Kf S „ K rs ) 
front and rear damper coefficient ( a/,, a r ) 

Front and rear tyre stiffness (K ft , , K rt ) 

Axle to axle distance (L) 

Polar moment of inertia of the sprung mass (J) 

The two-dimensional model of the car suspension system has four degrees of freedom, namely 

• Vertical motion of front unsprung mass( q{) 

• Vertical motion of sprung massf^) 

• Angular motion of sprung massO^) 

• Vertical motion of rear unsprung mass (^ 4 ) 
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Figure 2.9 Freebody diagram of front and rear unsprung masses and sprung mass 
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For these four degrees of freedom, the governing differential equations can be written as [7]: 
For vertical motion of front unsprung mass 

dt 1 U„ 

For vertical motion of sprung mass 


+ F 3 + F 4 + K;) 
dt 2 M s 

For angular motion ( pitching) of sprung mass 

d\ _ (F 4 + F s )L 2 - (F 2 + F 3 )L, 
Jt 2 J 

For the vertical motion of rear unsprung mass 

dt 2 M m 


where the Fi to F 6 are the forces due to relative deformation in the springs and velocities in 
damper ( Figure 2.9 ). These forces can be written with di , 62 , d 3 , d 4 , as relative deformations in 
front tyre, the front spring, the rear tyre and rear spring as, 


Fi-Kftd] 

F 2 =K fs d 2 


F 5 = <Xfd ,2 

F 4 -Krsd 4 

(2.5) 

F 5 = QLjdj 

F 6 =Krtds 
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Moreover these deformations can be written as follows 

di-qi-fi it) d 2 =q 2 +kqi-qi 

di=q4-f 2 {t) d4=q 2 -hq 2 -q4 ( 2 - 6 ) 

The time varying functions f\ (t ) and f 2 ( t ) are road irregularities on front and rear tyres as a 
function of time. As the car traverses a road, these roughness excitations at different axles are 
not independent. The rear wheel experiences nearly the same input profile as the front wheel, only 
delayed in time. The time delay is equal to the axle to axle distance divided by the speed of the 
travel. Solutions of equations (2.1) to (2.4) will give the bouncing and pitching dynamics of the 
car. 

The effect of road irregularities across the width of the car cannot be seen in the two-dimensional 
model. Also the above model cannot be used to get the dynamics of the car in rolling. Thus, we 
consider the three-dimensional model. 

2.3.2 Three-dimensional Model 

The need of complete analysis of the dynamics of the car necessitates the three-dimensional 
model for the suspension system. Clearly, the three-dimensional model should contain modeling of 
all four wheels of the car with the sprung mass resting on the two axles with primary suspension 
on all the four wheels. So, this model can be seen as a mass resting on four spring-damper pairs 
which are in turn, connected to a spring mass system. Such a model is shown in Figure 2.10. 

To define this model, in addition to unsprung masses, spring stiffnesses, damper coefficients, and 
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tyre stiffnesses at all the four wheels (front right, front left, rear right, rear left) we need the 
following considerations: 


• Polar moment of inertia of sprung mass in pitching ( J pitC h) 

• Polar moment of inertia of sprung mass in rolling ( J ro ii ) 

• Wheel base (axle to axle distance) ( L ) 

• Wheel track (distance between the right and left wheels) (W ) 

• Position of center of gravity of the sprung mass 

This model has eight degrees-of-freedoms. Governing differential equations can be written as 
follows: 


For vertical motion of front right unsprung mass 


d\ _ (F 2 +F 3 -F l ) 

dt M fur 


(2.7) 


For vertical motion of front left unsprung mass 

d\ = (F 5 + F 6 -F 4 ) 


For vertical motion of rear right unsprung mass 

d #3 _ (/g + ip — Fj) 

dt 2 ~ M mr 


( 2 . 8 ) 


(2.9) 


For vertical motion of rear left unsprung mass 


A 4 
dt 2 



( 2 . 10 ) 
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For vertical motion of sprung mass 

_ ~(^2 + ^ + F 6 + ^8 + + ffl + ^ 12 ) 

For pitching motion of sprung mass 


( 2 . 11 ) 


d 2 q 6 __ (F 2 + F 3 )A+(F 5 + F 6 )I 1 -(F 8 + F 9 )I 2 -(^ 1 + ^ 2 )L 2 

dt J pitch 


( 2 . 12 ) 


For rolling motion of sprung mass 


£& (F 2 + F,)L, - (F s + FJL, + (F t + F,)L, - (F u + F n )L t 

dt J w ll 


(2.13) 


The springs and dampers have no lateral resistance and as no lateral force on the model is 
considered. The yawing motion, the rotation about vertical axis has been ignored, in the above 
three-dimensional model. 

Hence, equations (2.7) to (2.13) describes the dynamics of the model in bouncing, pitching and 
rolling. By imparting different excitation functions ( fi(t') to fait)) to four tyres, any road excitation 
can be simulated. 

The forces and relative deformation in the system are given by 

F; = Kftr dl F 2~ Kfsr d.2 F 3= OL f sr d-2 

F 4 =Kfti di Fs—Kf S i d 4 F g= (Xf S i d 4 (2.14) 


20 



Chapter 2 


II 

(C 

F 8 — K rsr d() 

F g — ( %rsr ^6 


Fjo—Krti d 7 

Fll = K rs ld8 

F ;2 =: OC rsl 


di = qi 

d 2 = qs - q6 h - qi h~q\ 



d 3 = q 2 -flit) 

di = q 3 - q 6 li - q 7 U-q 3 


(2.15) 

ds = qi - flit) 

d(> - q$ +q6 h - qi h-q 3 



di= q± -fAit) 

d% = qs +q6 h + qi U~qA 




To calculate the dynamic response of the car, the steps are as follows: 

1. Derive the excitation functions fft) to / 4 (f) best representing the road conditions to be 
simulated. 

2. Calculate the relative displacements of all the points in the suspension system by using 
equation (2.15). 

3. Calculate the resistive forces in all the flexible elements of the suspension system using 
equation (2.15) in equation (2.14). 

4. These forces can be substituted in equations (2.7) to (2.13) to get the governing equations of 
the suspension system. 

5. These equations now can be solved for finding the dynamic response of the car to the 
excitation given by the road in form of excitation functions /i(r) to /4(f). 
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2.4 Summary 

The primary function of the suspension system is to keep the tyres on the road surface and to 
provide a smooth ride for driver and the passengers. The suspension is designed to absorb and 
dampen the vibration due to road roughness before they reach the passenger compartment. The 
road roughness profile can be represented as the sum of sine waves varying in amplitude, 
frequency and phase. The quality of the ride is a subjective perception. We are using the 
transmissibility of the vibrations as the index of the ride quality. To unprove open ride quality 
and/or handling properties of the vehicle, different types of the suspension systems starting from 
dependent front suspension and dependent rear suspension to strut type of front suspension and 
independent rear suspensions are in use with many variations among them. In this thesis, we are 
assuming a car with both front and rear suspension as of independent type. 

The quarter-car model of the suspension system does not give much information about the 
dynamics of the car. The two-dimensional model developed describes the motion of the car in 
bouncing and pitching but it fails in providing the dynamics of the car in rolling motion. This 
shortcoming of the two-dimensional model necessitated the development of three-dimensional 
model of the car. The three-dimensional model of the car has seven degrees of freedom and has 
21 parameters. Equations (2.7) to (2.13) describe the complete motion of three-dimensional 
model. With the help of this three-dimensional model, motion of sprung mass in bouncing, 
pitching and rolling can be analyzed. The road irregularities as a functions of time are used as 
excitation function for the model. 
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In the next chapter, we present a overview of genetic algorithm with a brief discussion of its 
applications in engineering problems. 
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Genetic Algorithms : An Overview 


Man has been always trying to mimic the nature. He makes use of principles of nature in order to 
make his living better. Genetic algorithm is one such example [12,14], Towards a step in 
mimicking the evolution, the effort in genetic algorithm has been to find solutions through 
computers . This chapter gives a bird’s eye view of the technique of genetic algorithm. 

3.1 Genetic Algorithm : Genetics and Algorithm 

Once there was a monkey and after few eras there developed a man. This is an example of 
evolution. Natural evolution involves perseverance and inheritance of good characteristics and 
simultaneous creation of new and better ones. On the most fundamental level all these are 
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achieved through evolution of thread-like structure called chromosomes. The perseverance comes 
through selection and creation comes through mixing and recombination of genetic material in the 
chromosomes. 

In a genetic algorithm (GA), all these are achieved through [12]: 

• Selection 

• Crossover, and 

• Mutation 

To perform these operations a GA need a population. The population has a number of members 
which are called individuals in GA terminology. Individuals contains variables coded like genes in 
chromosomes. In a simple GA, variables are coded in the form of strings containing alleles ( 
either 0 or 1). 

3.1.1 Selection 

On the population of individuals, usually the first operation is selection. Each individual, 
depending open its alleles have some decoded value and an associated fitness. On the lines of 
Darwin’s principle of survival of the fittest, this fitness acts as the criterion for the selection of an 
individual for taking part in crossover. The individual with higher fitness is selected more often 
then the one with lower fitness. Selected individuals are kept in a ‘mating pool’. In the mating 
pool, the individuals participate in creating new children points. 

For selection a number of schemes are developed, but the basic idea behind all of them is to pick 
above-average individuals for filling the mating pool with their multiple copies in some 
probabilistic manner. 
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One of the methods of selection is the stochastic remainder proportionate selection. In 
maximization problems, using this method, probability of selection of an individual is f, / Z f, , 
where f is the fitness value of the i-th individual and Zf, is the sum of the fitness of all the 
individuals in the population. Obviously an individual with large f, possesses higher chances of 
selection. This method is used for maximization problems. For minimization problems fitness 
function is to modified suitably [6,12]. 

Another most commonly used selection method is the tournament selection [6,12]. The advantage 
of this method is that it can be applied to both minimization and maximization problems without 
modifying the fitness function. In this methods “s” number of individuals are compared based on 
their fitness. For maximization problem the individual having highest fitness is said to be winner 
and selected to go to the mating pool while in minimization problems one with the lowest fitness 
value is said to be the winner and goes to the mating pool. The parameter, “s” is the number of 
individuals compared at a time and is called “tournament size”. Generally, the tournament size is 
kept as 2. 

3.1.2 Crossover 

“The remarkable ability of genetic algorithms to focus their attention on the most 
promising part of a solution space is a direct outcome of their ability to combine 
strings containing partial solutions” .[ Holland J.H. ‘Genetic Algorithms’ Scientific 
America July 1992] 

After creating the mating pool that is after selection, the second operator in GA is crossover. In 
nature, when sperm and ova fuse, mating chromosomes mix and exchange the genetic 
information. Similarly, in crossover operator, strings in the mating pool exchange information to 
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create new strings. Mating individuals are called parent points or parent individuals. In a GA two 
parent strings are picked up at random from the mating pool and across the length of the string a 
random cross site is chosen. To exchange the information, alleles on one site of the cross site are 
swept between the two individuals to create two new strings called “children points”. For 
example, if Si and S 2 two individuals are mating with a random cross site as shown 

si =110 0 11 

s 2 =10 110 0 

then after crossover the two children points s'i and s '2 will be 

s'i = 110 10 0 
s' 2 = 10 10 11 

The crossover is the operator which combines the genetic information of two solutions. The 
crossover is not done to all population members, instead it is done with a probability of p c ■ 

The above described crossover is called a single point crossover. There are many other cross over 
operators available, such as multi-point crossover, uniform crossover [12]. Many problem specific 
crossover operators also available [12]. 

3.1.3 Mutation 

After performing the crossover the new population is almost ready. In order to avoid premature 
convergence a mutation is used. As discussed in [12] “mutation is insurance against premature 
convergence”. This operator, “mutation” changes alleles ( 0 changing to 1 or vise versa ) 
occasionally with a probability p m . 

The mutation works as shown below: 
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1 0 0 Hi 0 1 — mutation— > 1 0 0 0 0 1 

The flow chart of a simple genetic algorithm is shown in Figure 3. 1. 

Mutation may provide a feature to the child string which none of its parents had. After doing 
mutation the new population is complete and a generation of GA is over. 


3.2 Applications Of Genetic Algorithms 

GA works through a coding of variables instead of variables themselves. So, GA works with a 
search space which is discretized, regardless of the function’s continuity. This increased the range 
of engineering problems that can be solved by GAs. With this the population-based approach of 
GAs enables the parallel processing of many solutions in different regions of the search space. 
Along with these advantages the fact that GAs requires only the fitness information from the 
objective function makes the GAs more robust then any other search and optimization algorithm. 
With these and many more advantages GAs had been successfully applied to a breathtaking range 
of problems, from design of mechanical components [8], structure designs [4] , metallurgical 
problems [18] and even to problems like criminal face recognition [2] . GAs are also used for 
solving optical filter design [5] and ambiguous shape modeling [1] problems. The list of 
application of GAs is growing very fast, GAs are getting more and more advanced with 
development of new ways of coding, new selection operators, new crossover operators and 
many other related developments. In this thesis, we use GAs to find optimal car suspension 
system for different objectives of a practical design. 
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3.3 Summary 

The search procedure based on the mechanics of the natural genetics and natural selection is the 
quintessence of genetic algorithms. The variables of the functions are coded in the form of a 
string. Strings get their fitness values from the underlying objective function. GA has three main 
operators, namely, selection, crossover and mutation. Selection is the process of forming the 
mating pool of successful individuals. In crossover, two individuals are chosen at random and 
substring on one side of a randomly selected cross site is swapped to create new children points. 
Mutation tries to introduce some extra information that is usually not possible by the crossover 
operator. Because of a number of advantages over other search and optimization algorithms, 
GAs have a wide variety of applications in many engineering fields. In this thesis we use GAs to 
find the optimal car suspension system. 
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Optimal Suspension Design Formulation 


In any system design problem, the objective is to first identify a number of salient parameters 
representing the problem and then to find a suitable set of values of these parameters. Different 
combinations of these parameters leads to different performances of the system. One may be better 
than others or some may be worse than other in respect to an objective. Thus there is a need to 
find that combination of system parameters which give the best performance of the system. One of 
the ways to achieve this is to pose the problem as an optimization problem. 

All possible parameters of a problem are not taken as design variables [7]. In an optimization 
problem salient parameters are first chosen as design variables. The other variables assumed as 
constant. Once the variables are identified, the next step is to chose an objective function. The 
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objective for which the optimization is to be carried out, should be very carefully chosen, because 
same problem will lead to different solutions for different objectives. The next step is to identify 
constraints which might arise to satisfy certain restrictions and limitations in behavior of the 
system. In this chapter, we discuss optimal design of the suspension system of a passenger car by 
identifying all components of a optimal design formulation, mentioned above. 

4.1 Two-dimensional Model 

The car suspension design can be achieved by considering two-dimensional and three-dimensional 
model of the suspension system. In a two-dimensional model the suspension system is modeled with 
only two wheels, one at front and other at rear. With two-dimensional model the vertical bounce 
motion and pitching motion of the sprung mass can be analyzed. As discussed in Chapter 3, equation 
(2.1) to (2.4) govern the dynamics of the four degree-of-freedom two-dimensional model of the 
suspension system. In the three-dimensional model, in addition to above stated motions of sprung 
mass, the rolling motion of the sprung mass can also be analyzed. We first discuss the two- 
dimensional model and then we will discuss the three-dimensional model. 

4.1.1 Mathematical Modeling 

Analysis of the equations (2.1) to (2.4) describing the dynamics of the two-dimensional model of a 
car suspension shows that they belong to the group of ordinary, nonlinear, second order differential 
equations. The solution of these equations in closed form is impossible to find. This is primarily 
because of the restrictions in spring and damper characteristics followed by the Tata Engineering and 
Locomotive Company (TELCO), Pune. For example, stiffness of rear spring and damping 
coefficients of dampers and at front and rear are not constant with respect to the deformation and 
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velocities which they are subjected to. This makes the coefficients of the equations variables and 
makes it difficult to find exact solution to differential equations.. Typical characteristics of the front 
damper is shown in Figure 4. la. Also, these elements have different resistance in bounce and recoil. 
To take this in account, the force equations in equation (2.5) are modified as follows: 

F 3 = C«/i +Yuv( a /3 -a /i)+Yi(l~YiwX a / 2 -cc fl )]d 2 -Yi(l-y 1N )(a /2 -oc n )a (4.1) 

F 4 =[K nl +y 2N (K„ 3 -K Bl )+y 2 {l-y 2lf )(K„ 2 -K nl )]d 4 ~Y 2 (.l~y 2N )(K n2 -K„ x )b (4.2) 

F 6 = [a rl +y 3N (a r3 -a rl )+y 3 (l-y 3N )(a r2 -a rl )]d 4 -y 3 (\-Y 3N )(a r2 -a rl )c (4.3) 

where y 2 is 1 for d 2 > a and zero elsewhere end yi N is one for d 2 < 0 and zero elsewhere. Similarly 
y 2 and Y 2 n are defined for d 4 and y 3 and y 3N are defined for d 4 . 

If the Z f , Y f , Zr , Y r , Zr S and Y rs are chosen such that 

F 3 = Z f d 2 - Y f a 

F 4 =Z rs d 4 -Y rs b 
F 6 = Z r d 4 - Y r c 

then to get the solution, the equations (2.1) to (2.4) can be transformed into a system of first- order 
differential equations. After the transformation they can be written as follows: 

Q = KmQ+C (4.4) 


Where 
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C = [K ft -m-aY f , aY f +bY rs +cY r , -b\J rs - cY r + alj f , -bY n - cY r + KJ 2 (t), 0, 0, 0, 


The equations can now be solved by using a numerical integrator such as Runga-Kutta method. The 
solution vector will contain eight elements corresponding to these eight rows of above matrix. Let 
the solution vector be { QuCh, Q3, Q4, Q5, (&, Q7, Qs } T , where Qi,Q 2 and Q 4 represent the vertical 
velocities of front unsprung mass, sprung mass and rear unsprung mass respectively. Q 3 represents 
the pitching velocity of the sprung mass. Q 5 , Q 6 and Q 8 represent the vertical displacements of front 
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unsprung mass, sprung mass and rear unsprung mass respectively. Q? represents the pitching 
amplitude of the sprung mass. 

4.1.2 Design Variables 

The two-dimensional model of the car suspension system given in Chapter 2 has eleven parameters 
in its description. But all of these are not chosen as design variables, some of them are really can not 
be the design variables. In order to make an optimization procedure converge in a reasonable time 
and keeping in view of their critical role in dynamic behavior of the system, following parameters are 
chosen as design variables: 

• Front spring stiffness (K fs ) 

• Front damper coefficient (ot f] ) 

• Rear spring stiffness ( K rs j ) 

• Rear damper coefficient (a rl ) 

In the previous section, it has been seen that, except front spring stiffness, all have different values in 
different ranges of displacements and velocities they are subjected to. Specifically, they can take 
three different values. To reduce the number of design variables, only one of the stiffness value and 
only one damping coefficient value for each damper is used as variable and the other two are 
assumed as linearly proportional to the variable. These two prepositional ratios will be fixed on the 
basis of the characteristics given by the ERC division of TELCO, Pune. 
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4.1.3 Objective Function 

In the simplest sense, the objective function is a mathematical relationship among design variables 
formulating the objective which is to be optimized. In this problem, a set of design variables gives a 
set of differential equations to solve and the solution of the equations gives the response of the 
system. From this response the objective is to be derived. The objective depends on the underlying 
road excitation. If the car is going over a bump, the shape of bump will act as road irregularity 
functions that are/; (t) and / 2 (f) in the formulation. The bump can be modeled as a half sine curve 
Then/// can be written as 

/ (f) = h sin(2itvt / /) 0<t< l/v 


Where 


h = height of the bump, 
v = velocity of the car. 

/ = breadth of the bump. 
t= time of travel. 

Rear wheel also experiences the same bump after time L/v , hence 


f 2 (; t ) = h sin(2ttv / lit - L / v)) L/v < t < (L+l)/v 


with L as the distance between rear and front axles. 

As the car goes over the bump, the sprung mass also goes up, but after coming back on level road 
sprung mass keeps on moving up and down depending on the underlying natural frequency of the 
car. This motion of the sprung mass gives a feeling of discomfort to the passenger. The amplitude of 
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this subsequent motion of the sprung mass should be minimized. In other words, the effect of bump 
transferred to sprung mass, that is the transmissibility of the displacement should be minimum. The 
transmissibility is used as the objective function for the two-dimensional model. Mathematically, the 
objective function is as follows: 

Minimize max(Q 6 ) lh (4.5) 

Since the sprung mass will follow the bump, during its motion over the bump, the transmissibility is 
defined after the respective wheel have passed the bump. 

4.1.4 Constraints 

Constraints are the criteria for a rejection or acceptance of a solution during the optimization 
procedure. The constraints thus play an important role in optimization procedure. 

The discomfort of passengers due to vertical motion of sprung mass does not end with its amplitude 
of the motion. The jerk, that is rate of change of vertical acceleration, experienced by the passenger 
also gives much discomfort to the passengers. So, for having comfortable ride, the jerk should be 
within prescribed limits. This constitutes one constraint for the car suspension design problem . It is 
a common practice in automobile industries to limit the jerk at 18m/s 3 [19]. Thus mathematically 

< 18 m/s 3 

or, 18-26 >0 (4.6) 

Complete NLP problem can now be stated as follows: 

Minimize max(Qe)//i 
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Subjected to 

18-&>0 
0< K& < 5.0 

0< K rsl , <5.0 

0<Ofi < 5.0 

0< OCri <5.0 

4.2 Three-dimensional model 

Three-dimensional model of the suspension system of a car enables us to include the rolling motion 
of the car. Effect of undulations of the road across its length and along its breadth both can be taken 
care of with the three-dimensional model. The seven degrees-of -freedom model (Figure 2.10), 
gives the description of the bouncing motion, pitching motion and rolling motion of the sprung mass. 
Motion of unsprung masses are also described by the model. Here we discuss this three-dimensional 
model of the suspension system. 

4.2.1 Mathematical model 

The three-dimensional model and its governing equations can be analyzed along similar lines to the 
above two-dimensional model. Differential equations are of the same nature and so are the flexible 
element characteristics. Mathematical complexity is increased in terms of number of equations and 
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number of terms in the equations. Due to nonlinear characteristics of the dampers and rear springs 
force equations in equation (2.15 ) should be modified as the following: 

^3 = + 7 in (tt fr 3 ~ a fri ) + 7 1 (1 ~ Yin - a frl )]d 2 - y x (1 - y lN )(a fr2 -a frl )a r (4.7) 

F 6 = + 7 2N(oc fl3 -a fll )+ y 2 ( 1 - y 2 N )(a )]d 4 - y 2 (1 - r 2 ^)(« n ~ « (4.8) 

■^8 = [F rsr i "I” Y 3N (K rsr2 — K rsrl ) + Y 3 (1 ~ YiN )(F-rsr2 ~ ^rsrl )]^6 — Y 3 (1 — YlN )(^ rsr2 ~ K rsr y )b r (4.9) 

F 9 =[CCrn + Y^GCm -«M) + r 4 ( 1 -r4N)(«^2 -Y ^~Y AN)(.Un 2 - a m)c r (4-10) 

F. i = [K rsll +Ysn (F rs i 3 — K rsll ) + Y s (1 ~ Y 5N )(F rs[2 — K rsll )]J 8 — Y 5 (1 — Ysn X-^ra/2 ~ F rsil )b t (4. 1 1) 

F n =[a rll +r 6 iv(a rf3 -<XHi) + Y 6 (X-Y 6 N)(<*ri 2 -Ye^-YeN^m -u r n)Ci (4.12) 

A^Hin WC CllOOSC Zfr , Zfl y Yf y > Yfl y Zrr r Yrr y Zrl , Yrl y Z rsr , Yrsr y Zr S i and Yrsi such that 

F 3 — Zf r di - Yf r a r F 5 = Zfi d. 4 - Yfi ai 

Fg= Zrsr d'6 “ Yfsr br F )/— Zrsl dg~ Y r sl bl 

F 9— Zrr dg - Yrr Cr F 12~~ Zrl d8~ Yrl Cl 

Transformation of equations (2.7 ) to (2.13 ) into a first order equation system leads to following 
matrix equation: 

Q = K U Q + C (4.14) 

Where 
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where 

Cz ~ Kfii ~ f%(t) — tyYfl 

C^K^-m-C^-bXsr 

C 4 = K nl -f A {t)-c l Y rl -b l Y rd 

Cs = a X/r 4 " a Xfi + b r Y rsr 4 " bi Y rd 4 - c r Y n + c l Y rl 

C 6 — —O-iYflly + biY rs 42 + cXX + C r Y r rll 4 " bfYrsrh. ~ t^Xfrk 

C 7 = a Xfll 4 + b t Y rsl l 4 + c l Y rl l A — c r Y rr l 3 — b r Y rsr h ~ a Xfrh 


These equations now can be solved by using Runga-Kutta method. The solution vector 
will contain fourteen elements corresponding to each row of above matrix. Let the solution vector be 
{ Qi >Ch , Q 3 , • • • •> Q 14 } t , Qi to Q 7 gives the velocities respective to {qi , q 2 ...qg} and Qg to Q w 
give the {qi ,q 2 , ,q 8 ) values. 

4.2.2 Design Variables 

Twenty-one different parameters are required to define the three-dimensional model of a car 
suspension (Figure 2.10 ). All of them are not taken as design variables. Springs and dampers of the 
suspension system are the elements used as design variables. 
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There are four springs and four dampers in the three-dimensional suspension system model. Two are 
at the front and two are at the rear. Spring and damper on front right side are considered to have 
exactly the same characteristics as the spring and damper on the front left side. Similarly the spring 
and damper at the rear left side are considered same as the spring and damper at the rear right side. 
Thus, in the 3-D case, the design variables are same as in the 2-D case. Explicitly stating, the design 
variables chosen for the problem are as follows: 

i) Front right spring stiffness(A/ jr ) 

ii) Front right damper coefficient a/ r ;) 

iii) Rear right spring stiffness(AT„ r ;) 

iv) Rear right damper coefficient a r w) 

For variables (ii), (iii) and (iv) values in different ranges will be taken in constant ratio based on that 
supplied by TELCO, Pune. 

4.2.3 Objective Function 

In the three-dimensional model, bouncing, pitching and rolling of the sprung mass can be analyzed. 
So, the problem can be formulated by taking any one of them on prime importance. Therefore three 
different objective functions are formulated to optimize design for each motion. The three objectives 
are as follows: 

• Minimization of the transmissibility 

• Minimization of amplitude in pitching 

• Minimization of amplitude in rolling 
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4.2.3.1 Minimization of Transmissibility 

For the purpose of the minimization of transmissibility, car’s motion over a speed breaker is 
simulated. The left and right tyres of the car moving over a speed breaker experiences the same road 
profile, with only the difference that rear wheel passes the same profile on the road a little late then 
the front tyre. In this situation, the forcing function f x (t) and f 2 (t) becomes the same. As the speed 
breaker is modeled as half sine curve, these functions can be written mathematically as follows: 

fi(t)=hsin(2Kvt/l ) 0< t < l/v 

m=m 

and also 

fi(t)=hsin(2nv/l(t-L/v )) L/v <t< ( L+l)/v 

m=m 

With these forcing functions, the response of the system can be found by solving equation (4.13), 
and objective function for the purpose can be written as 

Minimize max(Qi 2 )//i (4.14) 

Since the sprung mass will follow the bump, during its motion over the bump, the transmissibility is 
defined after the respective wheel have passed the bump. 


4.2.3.2 Minimization Of Amplitude In Pitching 

When a car traverses over a road having large undulations, high pitching of sprung mass may result. 
In such cases the objective of the suspension design should be to minimize the pitching amplitude. 
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To simulate such a road, a sinusoidal road with wave length equal to the wheel base is simulated. 
Road excitation functions for this type of road can be written as follows: 

fi(t)=hsin(Kvt/l ) 

m=m 

f 3 (t) =hsin( ( kv/1 ) (t-LUv) 

m=m 

Solving the equation (4.14) enables us to write the objective function as 

Minimize max(Qj 3 ) (4.15) 

4.2.3.3 Minimization Of Amplitude In Rolling 

If the right and left wheels moves over different road conditions, rolling motion will take place. Most 
severe rolling will be encountered when one wheel is going over a bump and other is going in a hole 
and vice versa. For getting the objective function, car’s travel over a twisting road in simulated. A 
twisting road in one which has waviness and left and right wheel have out of phase excitation. 
Excitation functions for such a road can be written as follows: 

fiit) = hsininvtll) 

m=-m 

fi (t)= hsin{(nvll)it-L/v)) 

Mt) = -MO 

Equation (4. 13 ) now can be solved with these forcing functions. The objective function then is 

Minimize max(Qu) (4.16) 
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4.2.4 Constraints 

Constraints for the optimization problem come from other related aspects which are not covered in 
the objective function formulation. In order to make the suspension design more realistic and 
practically applicable, a few other constraints are induced based on discussions with engineers of 
TELCO, Pune [19]. 

i) As discussed in the two-dimensional model, the jerk experienced by the passenger should be 
within prescribed limits. Mathematically, it reduces to the following: 

gi s 18-Q 12 >0 (4.17) 

II) Natural frequencies of the various subsystem affects the performance of the car. The bounce 
frequencies of front and rear suspension and that of pitching are most important. These frequencies 
can be calculated by using the Guest model [3] (Figure 4.2 ). Based on this model, the frequencies 
are given as follows: 

Front suspension frequency in bounce 


1 I 2K fi (a - x + j) 1 2 + 2K ~ (s-b- x) 2 (4.18) 

Ff 2;t V M(p 2 +s 2 ) 

Rear suspension frequency in bounce 


1 /2^(a-x-r) 2 +2g ra (r + Z> + x) 2 (4.19) 

2tt V M(p 2 + r 2 ) 
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Frequency in pitching 

(4.20) 

where 

x- distance between center of gravity and spring center 
r, s = distance of front and rear conjugate points from spring center 
a, b= distance of the spring center from front and rear 
p= radius of gyration of sprung mass 

These can be calculated by using following expressions ( Refer Figure 4.2) 



= — 

r 


Now using the expressions (4.18 ) to (4.20 ), constraints based on frequencies can be formulated. 
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i) Front suspension frequency in bounce should be less than the rear suspension frequency in 
bounce [10,17,19]. This is because, it will make pitch motion to die faster. 

Ff < F r 

or g 2 =F r - F f >0 (4.21) 

ii) The natural frequency of human body lies between 1.5 Hz to 2 Hz [11,17], So any of the natural 
frequency of any subsystem in a car should not be greater then 1.5 Hz. This fact can formulated by 
putting the largest of these frequencies to be less then 1.5 Hz. 

max( Ff , F r , F p ) < 1 .5 Hz 

or g 3 s 1.5-max(Ff,F r ,F p )>0 (4.22) 

iii) To avoid a very small natural frequency any of these frequencies has been kept to be more than 
0.8 Hz. This can be formulated as constraint by putting least frequency to be greater than 0.8 Hz. 

min( F f , F r , F p ) >0.8 

or g 4 s min( F f , F r , F p ) - 0.8 >0 (4.24) 

This completes the formulation of the constraints. Thus , there are four constraints to the problem. 

4.3 Summary 

The optimal design problem or any optimizations problem can be seen in three part, namely design 
variables, objective function and constraints. Every parameter associated with the design can not be 
taken as design variable. Few parameters which are most critical for the design are to be chosen as 
design variables. Formulation of the objective function requires the understanding of the problem so 
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as to choose an appropriate criterion as an objective function. The constraints but some sort of 
boundaries to the search space, as the acceptance of a feasible solution depends on them. 

Objective functions for both two-dimensional and three-dimensional models are formulated for 
solving the governing differential equations in a reduced form. Two-dimensional problem has one 
constraint, restricting the jerk experienced. The objective is to minimize the transmissibility of the 
road irregularity to the sprung mass. Design variables are the damper coefficients and spring stiffens 
at front and rear suspensions. 

Three-dimensional model has considered bouncing, pitching and rolling motion of the sprung mass. 
So, it can give more practically applicable results. The number of constraints for the three- 
dimensional problem are increased to four by putting constraints on the natural frequencies of the 
model. Three different objective functions are considered for these motions. The design variables are 
same as two-dimensional because the left and the right suspensions are usually the same. The NLP 
problem for the three-dimensional case can be written as follows: 

minimize Equation (4.15) ,(4.16) or (4.17) 


subjected to 


18-$*i2 >0 

1.5 - max( Ff , F r , F p ) >0 
F r - F f >0 


min( F f , F r , F p ) - 0.8 >0 
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0 < K fsr <3.0 
0 < K rsrl < 3.0 


0 <a frl < 3.0 


0<tVi <3.0 

In the next chapter, we solve the above two-dimensional and three-dimensional problems using a 
genetic algorithm. 
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Simulation Results 


In this chapter, results of the optimization carried out for the design of suspension 
system are presented. Optimization results for the two-dimensional model using traditional 
methods and using genetic algorithm are presented first. The proof-of-principle studies for the 
working of the code developed for the simulation and genetic algorithm are also presented. 
Thereafter, optimal design of the three-dimensional model is presented. Different motions of the 
sprung mass are plotted for the solutions found with different objective functions discussed in the 
previous chapter. The results are also compared with the design currently used at TELCO, Pune. 
Modeling of the a real road as a polyharmonic function is assumed and subsequent studies for 
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finding the optimal suspension design for maximum comfort as suggested in ISO 2631 are also 
discussed. 

5.1 Two-dimensional Model 

With the help of the two-dimensional model of the suspension system discussed in Chapter 2, the 
motion of the car over a bump (simulated with height 70 mm and breath 500 mm) is simulated 
and optimal suspension system design is achieved via solving the NLP problem discussed in 
Chapter 4. In this NLP problem, we have four design variables and one constraint. The objective 
is to achieve the minimum transmissibility of the road irregularities to the sprung mass. We solve 
the problem with some traditional optimization techniques first and thereafter with genetic 
algorithm. 

5.1.1 Traditional Methods 

A number of traditional methods are available for solving the optimal design problems. To find 
optimal design of the car suspension system, we used ‘constr’ routine of the ‘MAILAB 
OPTIMIZATION TOOLKIT’ [13]. This routine finds the constrained minimum of a function of 
many variables, starting at an initial estimate. Mathematically, the problem should be in the 
following form: 

Minimize f(x h ....x n ) 

Subjected to 

g } (Xi )<0 j=l to k 

To put our problem in this form the constraint on the jerk is rewritten as follow : 

#2 
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This routine uses a Sequential Quadratic Programming(SQP) method. In this method, a Quadratic 
Programming(QP) subproblem is solved at each iteration. An estimate of the Hessian of the 
Lagrangian is updated at each iteration using BFGS formula [7,13]. A line search is performed 
using a merit function. The QP is solved using an active set strategy. 

The governing differential equations of car dynamics are coded in MATLAB program min g 
language and a MATLAB routine ‘ode45’ is used to solve these equations to enable us to 
calculate the objective function and constraint values. The machine used for the purpose is a P.C. 
with PENTIUM processor running at 133MHz. Initially when the code was running with the 
default setting of accuracy in ‘constr’ , the machine took seven hours but was not able to solve 
the problem. In many occasions in these seven hours, it prompted the massage 

''Matrix is near singular or badly, scaled results may be inaccurate’ . 

Thereafter, when the option for the accuracy was set for a lower accuracy, the routine converged. 
In this way some results are found with different initial guesses. However, the results found by 
‘constr’ are very much dependent on the initial guess. Results with different initial points are 
shown in Table 5.1. 


Table 5.1 Results Obtained With SQP method of MATLAB 


Initial Solution 

Solution Found 

K fe 

Ofi 

Krsl 

Orl 

Trans. 

Const. 

K fe 

CXfi 

Krsl 

CXrl 

Trans, 

Const. 

1.0 


1.0 

1.0 

0.96 



0.23 

0.83 

4.99 

0.49 

-0.995 

1.5 

mm 

1.0 

0.75 

0.82 

-0.998 

1.85 

2.97 

0.1 

0.34 

0.57 

-0.998 

1.0 

1.0 

1.5 

0.9 

0.89 

-0.997 

3.76 

2.89 

3.10 

2.19 

0.48 

-0.996 


54 


































Chapter 5 


Since the SQP method converges to different solution when started from different initial solution, 
it can be concluded that either the function is highly non-linear and multimodal in nature or the 
algorithm is unstable and not efficient and suitable for solving this problem. Whatever may be the 
case, the algorithm also takes a considerable amount of time in converging to the solution. Thus, 
we decided to use GA, an efficient optimization method, which could perform well in nonlinear, 
and multimodal problem and does not take too many function evaluations to converge to a near- 
optimal solution. 

5.1.2 Genetic Algorithm (GA) 

Difficulties encountered m solving the NLP problem for the two-dimensional model with 
traditional optimization techniques turned us towards the genetic algorithm. Before starting the 
optimization procedure with GA, we perform a couple of proof-of-principle results. 

5.1.2.1 Proof-of-Principle Results 

Contour plots are extensively used in optimization studies to clearly show the feasible region and 
the optimal solution. By plotting any solution on the contour plot, one can easily verify the 
optimality of the solution. We follow this simple strategy to verify the algorithm and the code 
developed for the simulation of the car’s motion over the same bump, as used in the traditional 
case. It is easy to draw contour plots of the objective function on a two-dimensional plane. In our 
problem, there are four variables (K fe , On, Krsi, and a,i ). Since we can only vary two variables at 
a time to take advantage of counter plotting , there are total six combinations possible. We 
investigate two of these combinations. In both cases, the following procedure is performed: 

1. Only two of the four variables are varied and other two are kept fixed at certain values. 
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2. A contour plot of the objective function (transmissibility of the system in terms of vertical 
displacement) is first made by using an exhaustive search technique [7], 

3. The feasible region corresponding to the constraint on jerk is determined by using the 
exhaustive search technique and plotted on the contour plot of the objective function. 

4. From the feasible region on the contour plot the optimal solution is located by inspection. 

5. The optimization code is used for two variables only to find the optimal solution and the 
corresponding optimal objective function value. 

6. The optimal solution found using the optimization code is plotted on the counter plot to 
investigate the efficiency of the algorithm and working of the code. 

5.1.2.1.1 Varying Front Suspension Parameters 

In the first case, front spring stiffness ( K& ) and front damper coefficient (ocn)are varied and the 
rear suspension’s spring and damper are fixed at following values: 

K rs i- 2.86 kg/mm a r i = 1.01 kg-s/mm 

Figure 5.1 shows the contour plot of transmissibility with respect to different front suspension 
parameters. The feasible region on the plot satisfies the constraint for the jerk. The Figure 5.1 
shows that the solution space (feasible region) is disjointed. It has some infeasible region bounded 
by feasible region and vise versa. Multimodality of the feasible region comes from this 
disjointiveness and makes the problem difficult for traditional methods. The GA is used next to 
find the optimal solution of this two variable problem. The parameter used for the GA run are as 
follows: 

Lower and upper bounds on variable one ( Kf S ) 0.1 and 5.0 
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Lower and upper bounds on variable two ( On ) 
Population size 

Chromosomlength for variable one ( K fe ) 
Chromosomlength for variable two ( otn ) 
Maximum number of generations 
Crossover probability 
Mutation probability 
Tournament size for selection 


0.1 and 5.0 
: 30 
: 8 
: 8 
: 30 
: 0.9 
: 0.01 
: 2 


The optimal solution found by the GA is as follows: 

Kf S = 4.53 kg/mm, o tp = 1.72 kg-s/mm 

The solution found by the GA is plotted in Figure 5.1. The solution seems to be close to the true 
optimal depicted by the feasible region shown in Figure 5.1. Thus, it can be argued that GA has 
found a near-optimal solution. We vary the other two variables, next and repeat the above 
experiment. 

5.1.2.1.2 Varying Rear Suspension Parameters 

The variables for the front suspension’s spring and damper are fixed at following values: 

Kf s = 4.53 kg/mm OLp = 1.72 kg-s/mm 

The transmissibility value with respect to rear suspension parameters is shown as a contour plot 
in Figure 5.2. The feasible region on the plot satisfies the constraint for the jerk. The GA is used 
to find the optimal solution of this two-variable problem. The parameter used for the GA run are 
as follows: 
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Population size ; 30 

Chromosomlength for variable three ( K rs i ) : 8 

Chromosomlength for variable four ( 0 ^ ) : 8 

Maximum number of generations : 30 

Crossover probability : 0.9 

Mutation probability : 0.01 

Tournament size for selection : 2 


The optimal solution found by the GA is as follows: 

Kni =2.85 kg/mm, a r i = 1.02 kg-s/mm 

The optimal solution found by the GA and the optimal solution that can be visually inspected by 
from contour plot agree each other. Thus, it can be said that the GA has found a solution close to 
the optimal solution. 

The above results give us the confidence to use GA as an optimization technique for the car 
suspension problem having four variables. 

5.1.2.2 Varying Both Front and Rear Suspension Parameters 

In this ran now all the four variables are kept as decision variables. All other suspension 
parameters of the model are kept according to the data given by TELCO, Pune. These are as 
shown in Table 5.2. The car is assumed to ride over the same bump as in traditional case. 
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Table 5.2 Parameters for two-dimensional Model 


Front Tyre Stiffness 

15 kg/mm 

Front Unsprung Mass 

50 Kg 

Rear Tyre Stiffness 

17 kg/mm 

Rear Unsprung Mass 

115 Kg 

Sprung Mass 

730 Kg 

Axle To Axle Distance 

| 

2825 mm 

Distance of C.G. From Front 

1325 mm 

Moment of Inertia 

2.8944e+04 kg -mm 2 


The parameters used in this run are shown below. A population size of 30 and maximum of 30 
generations are found to be adequate for this problem. The string length for each variable is kept 
as 10, so there are 2 10 or 1024 values possible for each variable. This can give (2 10 ) 4 or 1.1*10 12 
different solutions for the problem. The crossover and mutation probabilities are kept in 
accordance with suggestions in GA literature [12]. 


Population size : 30 

Chromosomlength for variable one(K fe ) : 10 

Chromosomlength for variable two(otfi ) : 10 

Chromosomlength for variable three (K rs i ) : 10 

Chromosomlength for variable four((Xri ) : 10 
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Maximum number of generations 

: 30 

Crossover probability 

: 0.9 

Mutation probability 

: 0.01 

Tournament size for selection 

: 2 


GA is run on a HP-9000 machine, using ‘UNIX C’ compiler. Five different runs were performed 
with different random seed numbers. Figure 5.3, shows the progress of the population-best 
solution with number of generations. The figure shows that the best transmissibility reduces 
continuously with generation number. The optimal solution found has a transmissibility is equal to 
0.32 and the suspension parameters are as follows: 

Kf S = 4.53 kg/mm a,fj = 1.72 kg-s/mm K rs i = 2.86 kg/mm a r] = 1.01 kg-s/mm 

Time taken to find the solution in 90 seconds ( CPU time) .The constraint is not violated in this 
case. Since, with four variables contour plots can not be achieved, we show the superiority of the 
solution by comparing its performance on the dynamic response with that used by TELCO, Pune: 

Kf S = 1. 56 kg/mm a fI = 3.3 kg-s/mm K rs i = 1.45 kg/mm a r i = 1.00 kg-s/mm 

With these parameters the car motion is simulated for the first 10 seconds after the hitting of the 
bump by the front tyre. In Figure 5.4, the vertical motion of the sprung mass (bouncing motion) 
is plotted and in Figure 5.5, the pitching motion (angular displacement) of the sprung mass is 
plotted. 

From the Figure 5.4 and 5.5 it can be easily seen that the vertical displacement through out the 
10 second simulation (and hence the transmissibility ) is less with the parameter values resulted 
from optimization through GA. 
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Figure 5.3 Progress of GA for Minimization of Transmissibility (2-D model) 
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The value of displacement transmissibility (Eq. 4. ll)for the motion over a bump of 70 mm 
height and 500 mm spread is 0.82 with TELCO design. The same is reduced by about 62% to 
approximately 0.32, for the optimal suspension design. This reduction in transmissibility shows 
the effectiveness of the optimization procedure which is carried out with the help of GA. 

With these encouraging results, achieved through two-dimensional model optimization, we carry 
out the optimization, and present simulation results for the three-dimensional model. 


5.2 Three-dimensional model 

To study the effect of various road conditions on suspension design, we have used three 
different excitation functions and three different objectives. In all the cases, the following GA 


parameters are used: 

Number of runs : 5 

String length for variable K fe : 10 

String length for variable On : 10 

String length for variable K rs i : 10 

String length for variable 0Cri : 10 

Population size : 30 

Number of generations : 30 

Crossover probability : 0.8 

Mutation probability : 0.01 

Tournament size for selection : 2 


The best solution found in these runs is reported as the optimal solution. Simulation runs over 
different roads are presented next. 
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5.2.1 Over a Speed Breaker 

In this case the car is assumed to move over a speed-breaker. The speedbreaker is modeled as a 
half sine wave with following specifications: 

Height of the bump = 35 mm 
spread of the bump =100 mm 

The GA is then used for minimization of the transmissibility defined on vertical displacement (Eq. 
4.14 ), with constraints on jerk and frequencies ( Eq. 4.17, 4.21, 4.22, 4.23). The progress of GA 
with generations is shown in Figure 5.6. The best solution found in the initial population had a 
non-dimensionalized transmissibility equal to 0.48. After 30 generations, the optimal solution 
found is as follows: 

Kf S r= 1.54 kg/mm <X/ r i = 2.63 kg-slmm K rsr i = 1.29 kg! mm a rr i = 2.05 kg-slmm 
The corresponding value of transmissibility is 0.44. The bouncing motion of the sprung mass is 
plotted in Figure 5.7 and the pitching motion is plotted in Figure 5.8. The TELCO suspension 
design is an infeasible solution for this excitation function as the maximum value of the jerk is 
much more than the limit. The constraint violation by this solution is for jerk constraint, this give a 
jerk of 27000 mm/s 3 . Since, the existing TELCO solution is not a feasible solution, it will not be 
justifiable to compare our solution with TELCO design. 



Different GA Runs for Minimization of Gransmissibility 



Fisure 5.6 Progress of GA for Minimization of Transmissibility (3-D model) 
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5.2.2 Over a Sinusoidal Road 

As discussed in section 4.2.3. 2, a sinusoidal road is a road with the profile resemb ling a sine wave 
of large wavelength. The motion of the car moving with a velocity of 5 kmph on this type of the 
road is simulated. The road is defined using following specifications for the sine wave; 

Wave length of the wave = 2850 mm 

Maximum amplitude of the wave = 35 mm 

The GA is then used for the minimizing the maximum pitching amplitude (Eq. 4.15 ), with Eq. 
4.17, 4.21, 4.22, 4.23 as constraints on jerk and frequencies. The GA progressed towards the 
optimal solution is shown Figure 5.9. The optimal solution found after 30 generations is as 
follows: 

Kf„= 0.84 kg/mm a pi = 1.77 kg-s/mm K rsrJ = 0.83 kg/mm a rr! = 1.21 kg-s/mm 

The Figure 5.10 shows that GA has reduced the pitching transmissibility over a few generations. 
All three motions of the sprung mass, namely bouncing, pitching and rolling with suspension 
system parameters for this design are shown in Figures 5.11, 5.12 and 5.13, respectively 
(Presented after discussion of the subsection 5.2.4). 
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Figure 5.9 Progress of GA for Minimization of Maximum Pitching Amplitude 
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5.2.3 Over a Twisting Road 

Rolling motion of a car will be most severe when the car is moving over a twisting road. On a 
twisting road, right and left wheels moves out of phase on two sinusoidal road profile. Car’s 
speed is again taken as 5 kmph and motion on such roads of the following sine wave 
specifications is simulated: 

Wavelength = 2850 mm 

Maximum Amplitude = 35 mm 

With the objective of minimization of the rolling amplitude(Eq. 4.16), GA runs are performed. 
Constraints on jerk and frequencies are taken as in the previous cases. The best solution found 
(which gives the maximum rolling amplitude of 0.032132 radians) is as follows: 

Kf S r= 0.64 kg/mm a fr! = 1.02 kg-s/mm K rsr i = 0.43 kg/mm a rr i = 1-33 kg-s/mm 

Figure 5.10 shows the progress of GA with the generation while minimizing the maximum 
amplitude of rolling. Bouncing, pitching and rolling motions of the sprung mass with above 
values as the suspension parameters are shown in Figures 5.11, 5.12 and 5.13 respectively. 
(Presented after discussion on subsection 5.2 4) 
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5.2.4 Over A Combined Road 

Above road profile^ bump, sinusoidal and twisting) are simple and makes the optimization 
procedure easier to converge. Therefore, they are not totally practical . A road can be better 
approximated by a combination of these three basic road profiles. Thus, we find the optimal car 
suspension for all the three road profiles simultaneously. The car’s motion is simulated over each 
of the them and corresponding objective function value is calculated. These objective functions 
values are then combined to get an overall objective for the GA. GA runs are performed with the 


following parameters: 

Maximum number of runs : 5 

String length for variable K fs : 10 

String length for variable a// : 10 

String length for variable K rsI : 10 

String length for variable a rl : 10 

Population size : 30 

Maximum number of generations : 40 

Crossover probability : 0.8 

Mutation probability : 0.01 

Tournament size for selection : 2 


The best solution found in these runs gives the following values of the suspension parameters: 
Kfsr= 1.43 kg/mm a fr i = 2.46 kg-s/mm K m j = 1.55 kg/mm a rr ; = 1.91 kg -s/mm 
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For this solution, the different objective functions have following value: 

Transmissibility = 0.48 

Maximum pitching amplitude = 0.0253 

Maximum rolling amplitude = 0.03231 

Figure 5 1 1 shows the bouncing motion of the sprung mass for the four solution found above- 
optimized for the bouncing only, optimized for pitching only, optimized for rolling only, and 
optimized for all of them together.. The lowest value of transmissibility, as expected is for 
minimization of bouncing case. The combined solution gives slightly higher value for the 
transmissibility. The other two solutions have even higher vertical transmissibilities. 

Figure 5.12 depict the pitching motion of sprung mass for all the four solutions. Maximum 
pitching amplitude remains more or less the same. It is slightly higher for the solution obtained for 
the combined run. 

Figure 5.13 gives the rolling motion of the. sprung mass for all of the solutions found above. 
Maximum rolling amplitude is more or less same in all cases but it is slightly higher for the 
solution obtained for the combined run. 

The above study for various objectives and road profiles reveal the fact that if optimization is 
done only for minimization of the transmissibility of vertical displacement, other two objectives 
and road profiles need not to be optimized. The optimization of vertical displacement 
transmissibility inherently optimizes other cases. Thus, designers of car suspension may only 
optimize for vertical displacement. 

In next section we present, more practical approach to the suspension design. 
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5.2.5 Realistic Road and Comfort Ride 

In practice, the roads are neither of the types discussed earlier. The road undulations vary 
randomly m sizes and shapes. To simulate such a road, researchers have assumed a combination 
of sinusoidal road excitation having different frequency, phase and amplitude [9,16,20]. For this 
purpose, a polyharmonic excitation function is chosen which can be expressed as follows: 

X n 

/(*) = £ A(X)sm[2%^-t-$(X)], (5.1) 

A, 

where 

A(X)= Amplitude spectrum of road roughness 
X = wavelength 
t = time 

<j>(A,)= phase angle dependent on X 

The wavelength and amplitude has a linear relationship [9], which can be expressed as follows: 

A(X)=A 0 + Bo 

where Ao and Bo are the coefficient depending on the type of the road. Wavelengths in a ordinary 
road generally vary between 100 mm to 5000 mm. Hence, value of Xi in equation( 5.1) is taken as 
100 and Xn is taken as 5000. 

The phase angle is defined by the following expression: 

<j)(A.) = 2% (RND - 0.5) , 

where RND stands for a set of random number in the interval [0,1]. Now, the excitation function 
for the front wheels can be written as follows: 
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and for the rear wheels 


m=f 2 (t)=f(t) 


fi(t)=f 4 (t)=f(t-Uv) 

A profile of such a road created with above function is shown in Figure (5 .14). 

Passengers and driver of a car moving over this road will be subjected to vibrations varying in 
frequency and amplitude. The extent of vibration imparted to the passenger determines the 
comfort of the occupants of the car. Various methods of rating the severity of the exposure have 
been developed for specific applications. But, none of this methods are universally accepted [15]. 
In view of the complex factors determining the human response for vibration, international 
standards has been prepared to facilitate the evolution and comparison of the data and to give 
provisional guidance for suspension design. The standard ISO 263 1 defines and gives numerical 
values for the limits of exposure of vibration transmitted to the human body in frequency range of 
1 to SO Hz. This limits are given for the use according to the three generally recognizable criteria 
of preserving comfort, working efficiency, and safety or health. The limits set according to these 
criteria are named respectively in ISO 2631/1 the 'reduced comfort boundary’, ‘fatigue-decreased 
proficiency boundary’ and the ‘exposure limit’. In the design of passenger accommodations, the 
‘reduced comfort boundary is recommended as per ISO 2631/1. 


These limits are specified in terms of vibration frequency, acceleration magnitude, exposure time 
and direction of the vibration. Our simulation of the car motion gives only the displacement and 
velocity of the responses with time. Hence, the acceleration is calculated by central difference 
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method and is Fourici transformed to get the frequency dependency of the acceleration. Now, the 
acceleration can be plotted with frequency in accordance with ISO 2631. One such plot for 
TELCO’s design is shown in Figure 5.15. The plot shows that there axe some distinct frequency 
levels at which the amplitude of acceleration is severe. Passenger comfort level increases with the 
reducuon in these acceleration amplitudes at different frequencies. We exploit this fact to 
formulate the objective function for the maximization of the comfort. The procedure for 
calculating the objective function for a set of suspension parameters is as follows: 

1. For a suspension system, the three-dimensional dynamic response of the sprung mass is found 
for motion over the random road profile. 

2. Acceleration of the sprung mass is calculated using velocity of sprung mass history by central 
difference method. 

3. Fourier transformation of the acceleration data is done to find frequency- acceleration data. 

4. Peaks of the frequency-acceleration plot are recorded. 

5. Sum of the peaks of acceleration for different frequencies is taken as objective function value. 
With this objective function and constraints as in previous cases (jerk and frequencies), the GA is 
used to gel a optimal design (for minimum objective function value). The GA parameters similar 
to the previous cases are used to get the optimal solution. The optimal solution found after 30 
generations is as follows : 

K fsr = 1.4 5 kg/mm <X fr i = 0.14 kg-s/mm K rsr i = 1-28 kg/mm a rr i = 0.11 kg-s/mm 

It is surprising that the damping coefficients values in this solution are very close to their lower 
bounds. This perhaps has happened because we don’t have any velocity related constraint. 
Dampers are sensitive to velocity and since none of the constraints directly relate to velocity, the 


83 



Frequency Spec t: 


in 



5 Frequency (He) 

Fiaure 5. 15 Vertical Accleration of Sprung Mass- Frequency Spectrum 







Chapter 5 


optimization does not find the need of using a stiff damper. The frequency verses acceleration plot 
for this solution is shown in Figure 5.16. Which shows that the optimized suspension has an order 
of magnitude less acceleration than the TELCO design on the random road. Calculating comfort 
boundary from ISO 2631 chart, it is found that the car is comfortable for driving of 2.5 hours, in 
compare to only one hour in the case of existing TELCO design. This is a considerable 
improvement in the design of the car suspension. Table 5.3 shows an overall comparison of 
optimal suspensions obtained for different objectives and different road conditions. 

5.3 Summary 

Finding of all the optimization runs with different road excitations and different objective 
functions can be summarized in the following: 

1. Traditional optimization methods take about seven hours on PC with a PENTIUM processor 
and still do not found a stable solution. 

2. Optimal design of suspension found using GA, for 2-dimensional model gives 62% better 
performance than present TELCO design, in terms of vertical displacement transmissibility. 
Constraints on frequencies have some violation for this solution (table 5.3), but actually there 
were no such constraints on the design in 2-dimensional problem formulation. 

3. It was found that TELCO design does not satisfy the constraint on the jerk, in speedbreaker 
case with 3-dimensional model. Optimal solution found using GA is feasible and takes only 90 
seconds (CPU time) on a HP-9000 machine. 
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4. On a sinusoidal road, slightly better performance is achieved using GA, but the solution is 
closer to the optimal solution compared to the TELCO design, since the constraint values are 
close to zero. 

5. On a twisting road, considerable improvement in the rolling performance is achieved using 
GA. The third constraint on the frequencies has a value of -0.002, because the natural 
frequency of the front suspension in bouncing is 0.79799 Hz, which is very close to 0.8 Hz, 
and as the constraints on the lower and upper bound of the frequencies are relaxible, this very 
small constraint violation does not makes the solution infeasible. 

6. Optimal solution for the combined run shows a small violation of the constraint on upper 
bound on the frequencies because the natural frequency of the rear suspension in bouncing is 
1.515 Hz for the solution found. This values is very close to 1.5 Hz. Again, as the constraints 
on the lower and upper bound of the frequencies are relaxible, this small constraint violation 
does not makes the solution infeasible. 

7. On a random road, simulated by a polyharmonic function, the TELCO design fails to satisfy 
the jerk constraint. The maximum value of the jerk for TELCO design is approximately 
45000 mm/s 3 . The solution found using GA satisfies all the constraints and hence is a feasible 
solution. 

This thesis considered many objectives and road conditions and found optimal suspension. Since, 
the optimal solutions vary for different objectives, the designer may choose the one that is most 
suitable for application. Most importantly, this thesis has demonstrated that any such objective 
and road condition can be considered in the optimal formulation of the car suspension designed 
the optimal solution can be found out within a reasonable amount of time. 
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Chapter 6 

Conclusions And Extensions 


Based on this study we make some conclusions and then present some extensions. 

6.1 Conclusions 

Suspension systems of any vehicle is primarily used to isolate as far as possible the passenger 
compartment from the road undulations. Application of vibration theory helps us in understanding 
the vehicle dynamic response of a car to the road excitations. Transmissibility of the displacement 
has been considered as a criterion for the objective measurement of the vibrations. Two- 
dimensional as well as three-dimensional dynamic model of the car has been developed with front 
and rear both the suspensions as of independent type. 

In this thesis the work has been aimed towards finding the optimal suspension design. It has been 
observed that the optimal suspension design problem is too complex and multimodal to solve 
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using a traditional optimization technique borrowed for MATLAB OPTI MIZA TION 
TOOLBOX. The time taken by the machine was also considerably large. Genetic algorithm is then 
used to find the optimal solution. A number of different objectives have been used and a number 
of constraints are also checked as an acceptance criteria for the optimal solution. One constraint 
limits the value of the vertical jerk experienced by the passenger and three other constraints are 
related to the natural frequencies of the subsystems of the suspension system. Optimization 
problem formed in this way has four design variables corresponding to the front and rear 
suspension springs and dampers. 

It was observed that for minimum transmissibility, the springs of the suspension system tend to be 
more stiff. This increases the natural frequencies of the front and rear suspension. But as in two- 
dimensional, case the jerk was the only constraint, optimal solution did not take care of natural 
frequencies of the system. 

The three-dimensional model has been used next to find the optimal suspension system. Three 
different objectives, namely minimization of transmissibility, minimization of maximum pitching 
amplitude and minimization of maximum rolling amplitude have been used separately and the 
optimal solution is found for each of them. Thereafter, a combined objective function with a 
weighted average of above three objectives has been used to get the optimal design. The results 
suggest that minimization of transmissibility alone may be sufficient, because optimization of 
rolling or pitching amplitude results in more or less the same solutions. This happens due to the 
constraints on the frequencies. For rolling amplitude to be minimum, the spring stiffnesses will 
reduce to very low value, but lower bound on the frequencies keeps the stiffnesses above a certain 
limit and in turn keeps the maximum rolling amplitude on a particular level. On the other hand. 
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constraint of rear suspension frequency to be more then front suspension frequency takes care of 
pitching amplitude. 

In order to model a realistic road, a random road, modeled as a polyharmonic function has been 
used to excite the car motion. In an objective of maximizing ride comfort defined as per ISO 
2631 standard, GA has been able to find a much better solution compared to the TELCO design. 

In a nutshell, the following conclusions of the study can be made: 

1. Optimal car suspension design problem ( 2-D or 3-D) give rise to a feasible search space which 
is highly complex, multimodal, and disjointed. 

2. Popular traditional methods are not suitable to solve such problem. Although in some cases 
parameters can be adjusted to make the algorithm converge, the converged solution largely 
depends on the chosen initial solution. Moreover, the time taken to solve such problem is 
enormous. 

3. Genetic Algorithms (GAs) can solve such problem quite efficiently, because of their population 
approach and use of only function values. Moreover, the time taken to such problem is also 
reasonable. 

4. GAs have been used to solve the three-dimensional model of optimal suspension design 
problem for different road excitations and for three different objectives — vertical 
displacement, pitching, and rolling of sprung mass. GAs can repeatedly find the optimal 
solution. It is observed that the optimal solution for vertical displacement minimization for a 
car’s ride over a bump is also nearly optimal for other objectives and for other road models. 
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5. A polyharmonic road model has been built to simulate a realistic road. GAs have been used to 
find optimal car suspension for maximum ride comfort as per ISO 2631. The optimal solution 
found is better than that used by TELCO. 

Over all, it can be concluded that GAs are suitable to solve complex problem of optimal car 
suspension design. The success of GAs in this problem suggests their immediate application in 
other problems of engineering. Specifically, the following extensions can be imme diately looked 
at. 

6.2 Extensions 

Findings and investigations during the course of this thesis work open up a number of avenues for 
immediate studies in the optimal suspension design. We present some of them in the following: 

1. A more accurate and realistic model of the car for the suspension related studies can be 
formulated by assuming rigid body properties of the sprung mass. This will' help in finding 
more realistic responses of different elements in the suspension system. 

2. The model developed by us is for independent type of suspensions at front as well as at the 
rear. A model can also be developed on similar lines with independent front suspension and 
dependent rear suspension. This will be applicable to a large number of vehicles. 

3. Imposing a few constraints related to limits on the relative velocities experienced by the 
subsystems of the suspension will help in analyzing the system in accordance with ISO 
standards. 

4. Displacement of the unsprung mass at each tyre demands for the space, which is termed as 
suspension working space, may conflict with other features of the vehicle design. This also 
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puts some restrictions on the available solutions for the suspension design. Additional 
constraints satisfying these restrictions may find more practical solutions. 
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